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Abstract 

We consider the bremsstrahlung energy loss of high energy partons moving in the quark- 
gluon plasma, at weak coupling. We show that the rates for these processes receive 
large 0(g) corrections from classical (nonabelian) plasma physics effects, which are 
calculated. In the high-energy (deep LPM) regime these corrections can be absorbed 
in a change of the transverse momentum broadening coefficient g, which we give to 
the next-to-leading order. The correction is large even at relatively weak couplings 
as ~ 0.1, as is typically found for such effects, signaling difficulties with the perturba- 
tive expansion. Our approach is based on an effective "Euclideanization" property of 
classical physics near the light-cone, which allows an effective theory approach based 
on dimensional reduction and suggests new possibilities for the nonperturbative lattice 
study of these effects. 



1 Introduction 



The phenomenon of jet quenching, or suppression of high-p^ hadrons in A + A colhsions 
relative to expectations from scahng of binary p + p colhsions, has been the focus 
of much recent interest in RHIC physics |I] [2]. Its theoretical description ([3] and 
references therein) is based on the theory of jet evolution in thermalized media, whose 
uncertainties it is thus worthwhile to seek to reduce, or at least, quantify. This requires 
the calculation of higher-order effects, which we propose to do in this paper in the 
regime of weak coupling. 

As established by a large body of work on the thermodynamic pressure [i] [S] 
[6] , finite temperature perturbation theory meets with serious convergence difficulties. 
Unless the strong coupling as obeys as ~ 0.1, strict perturbation theory in powers of g 
is unreliable. Such a behavior seems generic: it is also observed for the next-to-leading 
order (NLO, 0{g)) corrections to thermal masses [7] [8] [9], as well as for the only 
transport coefficient presently known at NLO, heavy quark momentum diffusion [10] 
(whose behavior appears even worse). 

Following Braaten and Nieto j5j, who studied the thermodynamic pressure, these 
large perturbative corrections can be attributed to purely classical (nonabelian) plasma 
effects. They have shown this by first making use of the scale separation gT <C 27rT to 
integrate out the scale 27rT, leaving out a three-dimensional effective theory ("electric 
QCD", or EQCD) describing the scale rriD ~ gT as well as more infrared scales. The 
claim is then that contributions from the scale 27rT, as well as the parameters of the 
effective theory, enjoy well-behaved perturbative series |5] [11]; all large corrections are 
included in the effective theory. Furthermore, by treating this effective theory nonper- 
turbatively using various resummation schemes [9] [12] or the lattice [13], reasonable 
convergence can be obtained down to T ~ 3 — ST^. 

It is natural to expect large corrections from gfT-scale plasma effects in other quan- 
tities as well. Unfortunately, for real time quantities such as most transport coefficients 
and collision rates, a resummation program similar to that available in Euclidean space 
has yet to be fully developed and applied. This is because the real-time description of 
plasmas requires the hard thermal loop (HTL) theory [H] (which in essence is classical 
(nonabelian) plasma physics [15], also known as the Wong- Yang-Mills system [16]). 
which is arguably more complicated than its Euclidean counterpart EQCD. 

In this paper, we aim to point out progress which can be made for a specific class 
of "real-time" quantities: those which probe physics near the light cone. This includes 
the collision kernel C{q±) that is relevant for the transverse evolution of jets, whose 
crucial role in the theory of jet quenching will be reviewed below. 

To explain the idea, we first observe that the soft contribution to C{q±) (that 
arising from soft collisions with q± ~ gT) is described by soft classical fields that 
are being probed passively by the high-energy jet passing through them. These soft 



1 



classical fields are the fields surrounding the plasma particles. At this point we observe 
that field components moving collinearly with the jet are not particularly important 

— the standard calculation of collision rates fTf] (or see eq. fl22|) below) reveals that 
the contributing particles move with generic angles in the plasma frame, with even a 
suppression for the ones coUinear to the jet (due to the reduced center-of-mass energy) 

— which implies that the result must be insensitive to the precise value of the jet 
velocity f 1. The trick is then to think of f as f = 1+e — which, though unphysical, 
doesn't affect the answer — thus making the hard particle's trajectory space-like. This 
makes Euclidean techniques directly applicable including dimensional reduction, as will 
be explained below, thereby dramatically simplifying the calculation. 

In this paper we will thus (analytically) compute the full 0{g) corrections to the 
transverse collision kernel C{q±), describing the evolution of the transverse momentum 
of a fast particle. The second moment of that kernel gives the phenomenologically 
interesting momentum broadening coefficient q, which we also compute at NLO. 

This paper is organized as follows. In section |2] we summarize our results and explain 
their relevance to jet quenching; in particular we discuss the relevance of the parameter 
q. In section |3] we explain our computational strategy and formalism. Details of the 
calculation of C{q±) and of its (ultraviolet-regulated) second moment q are given in 
sections m and m respectively. In section[6]we derive, at NLO, the relation between the 
collision kernel C{q±) for momentum broadening and that for jet evolution — which 
turns out to be identical to the leading-order relation — and we discuss certain operator 
ordering issues which could enter higher-order treatments. Finally, in Appendix |A] 
we relate our approach to a slight generalization of sum rules previously found by 
Aurenche, Gelis and Zaraket [T8] . 

Alternative estimates of q and of jet evolution, based on gauge-string duality (see 
for instance [HI [20] [21] |22] |23] |2l] |25] ) , will not be discussed in this paper. 



2 Results 



2.1 Collision kernel 



The main result of this paper is the full next-to-leading order {0{g)) (analytic) expres- 
sion (l20l) for the two-body collision kernel C{q±), defined as: 

C{q±), (1) 



describing the evolution of the transverse momentum of a hard particle (with E T). 

The 0{g) corrections to C{q±), which being due to gT-scale physics only arise for 
q± ~ gT -C T, are illustrated in fig. [H Both the LO and NLO kernels C{q±) are 
proportional to the (quadratic) Casimir of the gauge group representation of the jet. 
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Figure 1: LO and NLO collision kernels C{q±) = {2TTYdT / (Pq^ for a fast quark in QCD 
(with iVj = 3), for = 0.1 and = 0.3. For gluons the curves are to be multiplied 
by a (Casimir) factor 9/4. 

The "leading order curves" is based on the full (unscreened) expression f l22l) at hard 
momenta, multiplied by q]_/ [q^+rn]^) to make it merge smoothly with the analytic 
result (fTOj) at low momenta, following the prescription given in [26]. The "next-to- 
leading order" curves use the leading order curves plus C(g_L)^^'"°'* given in (l20l) . 

The NLO correction is already quite large for = 0.1, giving nearly a factor 
of 2 around q± ~ T. As discussed in the Introduction, this is consistent with the 
behavior observed for 0{g) effects in other quantities. At as = 0.3, a typical value 
used in comparisons with RHIC data (see e.g. [27]), it is clear that the strength of 
the correction has grown out of control, meaning that (presently unknown) yet higher- 
order corrections are most certainly also important (though our results suggest that 
the value of needed to fit the data might be significantly smaller than the estimate 
of [27]). 

An interesting by-product of the approach used in this paper is that it extends 
naturally to higher orders: it makes perfect sense to evaluate the gauge-invariant Wil- 
son loop ([9]) nonperturbatively within the Euclidean three-dimensional EQCD theory, 
for instance using the lattice. Although this may not include all 0{g'^) corrections to 
C{q_[_) (contributions from the hard scale 2'kT will be missed), by analogy with the 
works on the pressure discussed in the Introduction, these missing contributions can 
be expected to be numerically suppressecG. We leave to future work the study of this 
interesting possibility. 

^Their description could turn out be very complicated, though, because jet evolution at O(g^) 
should contain, among other things, the analog of the NLO vacuum DGLAP splitting amplitudes in 
the presence of the LPM effect (described below). Also, various effects involving the scale evolution 
of the medium constituents and coupling constant evolution should arise. 
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2.2 Application to Jet Evolution 



The dominant energy loss mechanism of high energy particles (at weak coupling) is 
bremsstrahlung (including quark-antiquark pair production), triggered by soft colli- 
sions against plasma constituents. The theoretical description of these processes, at 
the leading order in the coupling, is well-established [28] [29] [30] . Their duration tform 
depends on the energy of the participants, and can interpolate between the Bethe- 
Heitler (single scattering) regime tform ~ E/q\ ~ E/mj^ at energies E < T, and the 
Landau-Pomeranchuk-Migdal (LPM) [31] (multiple-scattering) regime at high ener- 
gies E ^ T, with tform ~ a/ E/q, in which destructive interference between different 
collisions plays a significant role. 

In all of these regimes, however, the description factors into a "hard" coUinear split- 
ting vertex (Dokshitzer-Gribov-Lipatov-Altarelli-Parisi, DGLAP vertex [33]), times an 
amplitude (wavefunction in the transverse plane) which describes the in-medium evo- 
lution of the vertex. The latter accounts for the collisions which trigger, and occur 
during, the splitting process [28] [29] [30]. The DGLAP vertices themselves only in- 
volves hard scale physics (in essence, they are Clebsch-Gordon coefficients between 
states of different helicities) and thus cannot receive 0{g) corrections; the NLO ef- 
fects, which come from soft classical fields with p ~ gT, are included in their dressing 
amplitude. 

In section [6] we discuss these amplitudes at NLO and show that the relevant (three- 
body) collision kernel factors as a sum of two-body kernels C{q±), exactly like the LO 
one does [28] [29] [SOj |32]. As a consequence, our results can be used to give a full NLO 
treatment of radiative jet energy loss; one must simply include the NLO shift (!20|) to 
the two-body kernel C(g_L) which serves as an input to these calculation^. 

2.3 Momentum broadening coefficient {q) 

When the effects of a large number of small collisions are added together, it is natural 
to replace them by an effective diffusive process. The diffusion coefficient relevant for 
transverse momentum broadening, g, is defined as the second moment of the collision 
kernel ([I]): 



The ultraviolet cutoff |g_L| < gmax is needed to deal with the weak power-law falloff 
C{qA.) ~ 9^T^ at large g_L, which leads to a logarithmic dependence of q on gmax- 
This is a leading order logarithm; below we shall comment on the value of the cutoff 
9max- Using our NLO kernel (l20l) we can calculate the expansion of q up to terms of 

^ For instance, one would simply modify "C((7±)" in [32], which is actually equal to C{qi_) / {g^CsT) 
in our conventions. 




(2) 
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order g"^ §: 

q ^ Ca 
g^CsT^ Qn 

NfTf 



+ 



log ( — ) + ® log - 0.068854926766592 . . . 

KmoJ C(2) V T / 

log (—) + ^TT^log i^) - 0.072856349715786 
yniDj 2C(2) V T / 
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with ^(NLO) ^ _|_ (^371-2 + 10 -4 log 2) ~ 2.1985 a constant calculated in section O 
characterizing the NLO correction to q, and mj^ = g'^T'^{Nc + NfTf)/3 the leading- 
order Debye mass. In QCD with Nf = 3 flavors of fundamental quarks, Ca = 3 and 
NfTf = 1.5. For a discussion of the leading order result and logarithms, we refer the 
reader to P. Arnold's work [26j, from which the high-precision numbers were taken. 

The series ([3]) is meant to represent the expansion in g of the area under the curve 
of plots such as flg. [H For = 0.1 and a quark {Cg = |) the area under the leading 
order curve in the flgure (up to gmax = 4T) would yield q^'~' ~ 2.60 (T/GeV)^ GeV^/fm 
whereas the flrst two lines of © give q^'^'^^ ^ 2.08 (T/GeV)^ GeVVfm. The NLO 
shift is Aq ?a 2.22 (T/GeV)^ GeV^/fm from the flgure, about a factor of two effect, 
and Ag*^ ^ 5.26 (T/GeV)^ GeVVfm according to (El). Thus the third line of (ED itself 
suffers from sizeable truncation errors compared to our full NLO result fl20|) . We would 
like to stress, however, that ( !20l) . and flg. [T], is not merely a simple truncation error 
from a lower-order contribution but represents a genuine NLO effects. 

As discussed in the preceding subsection, it would be premature to attempt com- 
parison of our q results with experimental data, since it is clear that (yet unknown) 
higher-order corrections should also be important at physically relevant couplings. It 
is also worth noting that different approximation schemes taking q as input (this ex- 
cludes the AMY scheme [301 ES], which uses the full C{q±)) when fltted to RHIC data, 
tend to disagree rather signiflcantly on its preferred value [3]; since a critical analysis 
of these approximations lies beyond our scope, this simply means it is not completely 
clear which experimentally-extracted value of q we should comparing with. 

It seems appropriate here to recall some subtleties associated with the phenomeno- 
logical parameter q, which do not arise if one instead works with the full collision kernel 
C{q±). First, the value of the cutoff gmax to be used in (El) is process-dependent: since 
the g_L>gmax tail of C{q±) describes collisions occurring on a flnite rat(^ ^(q±>gn^a.^) ~ 
g'^T^ / q'^s.x: weighting them with in (1231) ceases to make sense for ^^q^^q^^^) ~ ^jet, 
with tjet the jet's lifetime, to be replaced with a formation time tform for bremsstrahlung 
pairs in the context of energy loss calculations. Therefore, parametrically, one should 



■^Version 1 of the present paper contained a transcription error which is now corrected in ([3]) . The 
argument of the logarithm (incorrectly) read q* as opposed to mjj obtained from (j23p and (j24l) . 
"'l am indebted to G. D. Moore for this discussion. 
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setgmax ~ 
a/E/^ so g, 




:orm 



Second, the presence of the ultraviolet tail implies that collisions having ~ 
(Jmax (with gmax the physical cutoff as determined above), which are intrinsically non- 
diffusive, already contribute at the next-to-leading logarithm order to bremsstrahlung 
rates. That is, they contribute at 0{1) compared to the log-enhanced diffusive con- 
tribution ~ loggmax/"^D comiug from mo -C 5'± <^ Q'max- Therefore, approximations 
based on diffusive physics lead, at best, to expansions in inverse logarithms of the 
energy. Such expansions were systematically studied in [31], with the conclusion that 
formulae to next-to-leading logarithm can be trusted at least when E^et ~ lOT, with the 
inclusion of the subleading term (e.g. the constant under the logarithm) being manda- 
tory. None of the presently applied approaches which rely on g as a phenomenological 
parameter presently includes such subleading logarithms (see [31 [35] and references 
therein for an overview of these approaches) . 

3 Strategy: Space-like Correlators and EQCD 

In this section, we relate certain correlators at space-time separation (more precisely, 
correlators supported on space- like, and light-like, hyperplanes of the type x° = vx^, 
V < 1), to Euclidean-signature correlators. We will then apply to them the formalism 
of dimensional reduction. 

3.1 Space-like correlators 

Field operators at space-like separated points (anti) commute with each other and their 
correlator does not depend on the operator ordering. For two-point functions at van- 



G>(t = 0,x) ^ 7^r^Tre-^^0.(x)O,(O)=T5^ /-^e^P-G^(u;.,p), (4) 



with the sum running over the Matsubara frequencies Un = 2iTinT, /3 = 1/T, and with 
Gfj the Euclidean correlator of some operators Oij taken here to be bosonic. 

In Lorentz-covariant theories, (jll) can be extended immediately to any correlator 
which is equal-time in a suitable boosted frame. Specifically, under 8b ^-cLXlS boost with 
velocity v, the thermal density matrix transforms to: 



ishing time separation, a well-known Euclidean representation holds 




n 




(5) 



^ This should be contrasted with the often-used kinematic cutoff 



max 



ET. 
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the primed quantities referring to quantities in the boosted frame; 7 = ^^^^-p- . The 
identification of H' and P'^ as the generators of time and space translation shows 
periodic identification x'^ = x'^ + i7(/3, —vP^^x) for the geometry associated to ([5]), 
with associated quantization condition on the "Matsubara frequencies" p'^ + vp'^ = 
27rinT/j. The spatial momentum p'^ must be kept real: it serves as a label for the 
physical states living on the x"^ = hyperplane. Thus only the frequency is complex. 
This determines the extension of (jl]) to equal-time two-point functions in the boosted 
frame: 

G>.(x'° = 0, x') = ^ 1 P'), P'n' = -vp" + 2mn^ . (6) 

It will be convenient to boost this formula back to the plasma rest frame, and to 
write it for general space-time arguments (w = ^): 



p° = 2mnT, pl=p^ + 2mnT—. (7) 







Eq. ([7]) is the main result of this section. It differs from (jl]) only due to the imag- 
inary part of p^, which, in fact, is required for the convergence of the sum over n: it 
prevents the Fourier exponential to contain exponentially growing terms as opposed 
to pure phases. Eq. ([7]) extends in a straightforward way to any higher-point correla- 
tor supported on = v)-type hyperplanes: one gets a summation-integration J2n Ip^ 
with the Pn as in (JTj), for all external legs, subject to the usual restriction of momentum 
conservation (and thus of "n conservation"), as for equal-time higher-point correlators 
|36j . The momenta running in loops must also be "twisted" like those in ([7]), e.g. 
Imp^ = vlmp^, to reflect the boosted- frame origin of the formulaic. This ensures that 
the imaginary part of every momentum is time-like, which is the natural domain of 
Euclidean physics. 

We will be interested in the amplitudes of ultrarelativistic dipoles moving with 
velocity v = 1, e.g. = x^ (note v = 1/v in general). Our derivation of ([7]) might 
seem compromised, since an "infinite" boost with velocity v = 1 obviously doesn't 
exist. However, a more careful look at the argument reveals that the boost plays no 
important role: after all we un-did it in the end. All that is really important, is that 
we can imagine quantizing the system along hyperplanes parallel to and express 
the thermal density matrix within these hyperplanes. Since it is certainly possible 
to quantize a system along light fronts, the result ([7]) must hold for x^ = x^. An 
alternative derivation of this, based on sum rules, is also given in App. lAl 



Note that this gives an impUcit dependence on the velocity v to GfAp'^,p^,p±^) appearing in ([7|). 
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The attentive reader might complain that setting = in ([7]) corresponds to 
taking a w \ 1 hmit, whereas the physically relevant regime v 1 lies beyond the 
reach of ([7]). Are we claiming that these limits are equivalent in general? Nc0. Our 
claim, explained in the introduction, is merely that for classical plasma physics effects 
they are equivalent — to the extent that the observable of interest is only passively 
probing a soft classical background, only a phase-space suppressed fraction of which 
propagates collinearly with the jet, this seems to be rather robust. It is unclear whether 
this will remain true when quantum effects are included (which will enter at ~ g'^), 
though, because of coUinear components present in the jet's own wavefunction. 

3.2 Dimensional reduction 

Naturally, the contribution from soft physics (momenta ~ gT) to sums like ([7]) is 
expected to be dominated by the n = mode. We will thus begin by "integrating out" 
the modes with n 0. 

First we claim that loop diagrams for which all external momenta have tt, = are 
equal to the standard ones. These two sets of diagrams have p° = and pz real, 
and could only differ due to the "twisted" Matsubara momenta ([7]) which circulate 
in the former; our claim is that this does not affect their value. The reason is that 
the imaginary part of every momentum P in such loops is time-like with a real part 
obeying Rep° = 0, ensuring that ReP^ is positive-definite. The imaginary part of 
the integration contours can thus be deformed from Imp^ = vlmp^ to Imp^ = 0, 
without crossing any pole. 

In particular, the modes with n = are described by precisely the standard "elec- 
tric QCD" (EQCD) three-dimensional effective theory [371 [38]. EQCD is pure three- 
dimensional Yang-Mills with coupling constant g^ = g^T coupled to a massive adjoint 
scalar Aq of mass mo- It is an effective theory for the gT scale (the Euclidean version 
of the hard thermal loop theory [H]), in which the loop expansion proceeds in powers 
of g'^T/niD ~ g. Its parameters do not receive 0{g) corrections. 

^ In strongly coupled theories accessible to gauge-string duality, these two limits are known to be 
physically distinct. A calculation of q for a physical massive quark moving with v < 1 (in the sense 
of its momentum broadening coefficient) by Teaney and Casalderrey-Solana [21], and by Gubser |22) . 
found a divergence q ^ (1 — w^)^^/'*\/AT^ as u 1. This calculation is valid for energies E < M^/XT'^, 
beyond which the coherence time of the force acting on the quark becomes of order the time scale of 
Langevin dynamics 22 ; a time-independent description is then impossible. This suggests that q for 
V < 1 should depend on a cutoff time scale, as is the case at weak coupling. 

On the other hand, the v\l limit has been studied by Rajagopal, Liu and Wiedemann [T9 l I20 j . 
by embedding Euclidean worldsheets into AdSs space, and no divergences were met in this limit. It 
is thus qualitatively quite distinct. 



8 



->x^=x3 



Xx 



Figure 2: Wilson loop representation of the dipole amplitude. 



The propagators of EQCD are: 

= G^q) = ^ - (8) 

+ mf) q^ 

(We use the tilde to denote that these are three-dimensional propagators.) The minus 
sign in front of the propagator reflects the fact that we will couple it to Minkowski- 
space Wilson lines: we have not performed a Wick rotation. 

In addition to its interaction with the n = modes, we must also include the 
direct coupling of the operator of interest to the n ^ modes. Physically, and as 
shown in Appendix [Xj a contribution from these modes would correspond, in the 
real-time formalism, to a failure of the soft approximation nsip^) ~ T/p^. Such a 
failure would signal a contribution from the T region in Minkowski space, which 
would necessary be signaled by ultraviolet divergences in the soft approximation, since 
this approximation correctly describes the intermediate region gT <^ p° <^ T and any 
contribution from the scale T should leave an imprint on this region. Thus, provided we 
do not find ultraviolet divergences from the n = contribution alone (which computes 
exactly the soft approximation, see Appendix this argument shows that we can 
safely ignore the direct coupling to the n ^ modes. This will turn out to be our case. 



4 The calculation 

In this section we express the collision kernel C{q±) as a correlator supported on = x° 
trajectories and evaluate it using EQCD. We only give details in the Feynman gauge 
^ = 0, though we have explicitly verified the ^-independence of our (gauge- invariant) 
collision kernel, check on the calculation. 



4.1 Operator definition of C{q±) and leading order result 

The evolution of the transverse momentum of a high-energy particle can be described 
by looking at its density matrix, as described in detail in [21]. For classical effects. 
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however (and even more so because we are taking the velocity to be f = 1+e), we can 
neglect operator ordering issues and replace the evolution of the density matrix by that 
of a dipole; we will come back to operator ordering issues in section I6.2[ 

High energy dipoles {E ^ mr)) propagate eikonally in the soft classical background. 
The collision kernel describing the evolution of its transverse momentum can thus be 
recovered from the Fourier transform of the (long-time limit of the) dipole propagation 
amplitude W 



^ C{q^) = [ d^x±e'P^-''^C{x±). (9) 



C{q±) is short for {2'7i)'^dT/d'^q±, as in ([T]). The dipole amplitude W(t,x±) is given 
by the trace of a long, thin rectangular Wilson loop stretching along the light-cone 
coordinate x"*", with a small transverse extension x± (see fig. [2]). 

The naive dimensional reduction of the Wilson loop iQ yields a Wilson loop stretch- 
ing along the 2;-axis of the three-dimensional EQCD theory. It couples to the linear 
combination = {A^ + Aq) of the EQCD fields, reflecting its ultrarelativistic origin. 
This "naive" dimensional reduction corresponds to keeping only the direct coupling to 
the n = modes. As explained in subsection 13. 2| this will be justified provided we do 
not find ultraviolet divergences. 

At the lowest order in perturbation theory, only the single-gluon exchange diagram 
((a) of fig. [3]) contributes. 



C{q^) = g'TCs dz d'x±e'P^-''^G+4z,x±) 

J — oo J 

= g^TC,G++{q, = 0, gx) = g'TC, (\ - ] , (10) 



^ql ql+mlJ ' 

where we have used (|8]). The compact form ffTOl) was first obtained by means of sum 
rules by Aurenche, Gelis and Zaraket [18], which we show in Appendix |A] are equivalent 
to our approach. 

4.2 Diagram (b) 

At the next-to- leading order (one- loop), self-energy insertions to the single-gluon di- 
agram (b) contribute (we will often write "gj," for a three-vector with = 0, which 
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Figure 3: Tree and one-loop diagrams contributing to C{q±). 



should cause no confusion; "J^" is short for J ■^:;^)- 
C{q^)ii»/9'TC, 

9'TCa 



Each bracket includes the contributions of one fish and one tadpole diagram, while the 
last one also includes the ghost loop. 

The (linear) ultraviolet divergences in ( ITT]) are to be canceled by matching counter- 
terms that can be unambiguously calculated within the framework of dimensional re- 
duction [371 EH]- They merely represent the (hard thermal loop) coupling of the n 7^ 
gluons to the soft n = ones, e.g. the gluon contribution to the mass squared 
m|). The fact that the direct coupling to exchange gluons with g° = 7^ does not 
contribute to the divergences can also be checked explicitly, from the convergence, with 
respect to g^, of the real-time integral ( l22l) (this justifies making the soft approximation 
on g°). Thus the divergences in ( ITTj) do not signal the presence of "new contributions" 
beyond the EQCD effective theory, as discussed in section 13.21 

Employing dimensional regularization, the divergences simply go awa}|§ and the 
counter-terms are zero to 0{g) [38J. This way we obtain (all our arctangents run from 
to 7r/2): 

Ciq±)ib) _ -m^-2^^tan j j ^^-_-_tan (^^j 

q^T^CsCa 4n{ql+miy 32gi + Svrgi ^ ' 

^ The dimensionally-regulated integrals ([TT|) have poles in dimensions 2 and 4 but are finite and 
unambiguous in dimension 3. 



(5noo(g. 



5W%qA 



{ql+miy ql 
(2g± - 



2pl 



3pl + 2ql + 

p2(g_L-p)2 



Pz 



p^{q±-pY, 
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4.3 Diagram (c) 

Diagram (c) plus its permutation contribute: 



)(c) 



g^T^CsCA 



1 



P Lgi(pH"^D)(('?±-p)H"^|)) iql+ml,)p'^{{q^-p)'^+m 

1 



-tan-^ ( ^ 



tan 



+ 



2nqj 



(13) 



(14) 



In the Feynman gauge, there is no contribution involving only transverse gauge fields: 
such a contribution would involve the (trivial) zzz vertex. Eq. flT3|) is manifestly con- 
vergent. 



4.4 Diagrams (d)-(g) 

Our calculation is based on a quasiparticle expansion, e.g. we simply set on-shell the 
external legs of scattering diagrams. The relevant expansion parameter is g, e.g. the 
ratio of the scattering width ~ g'^T to the scattering's natural frequency scale m^,. 
Thus in evaluating the external state corrections (d) we need only keep those effects 
which are not suppressed by the smallness of the width. A narrow resonance being 
described by just its position and the total area under it, this means that diagram (d), at 
0{g), produces only mass-shell corrections and wave-function renormalization factors. 
The (here imaginary) "mass-shell" corrections have no effects: they are identical for 
the initial and final states, so the "energy" (read z-momentum) transfer is zero in any 
case. The wave-function renormalization contribution is given by an energy derivative 
of the eikonal self-energy, and (e) is unambiguous, yielding respectively (including all 
diagrams of similar topology): 



C{q±){d) 
Ciq±)ie) 



2g'T'CtG^^{q^) G^M 



d 1 



p dp^ p^ - ie 

1 „ . ^ , . f ^ , . 1 



2/T^C,(C, - -Ca)G+M / G+4p) 



(15) 



The sum of (d) and (e) is proportional to Ga and identically vanishes in the abelian 
theory (Ga = 0), as required by abelian exponentiatioii^. 



^ Abelian Wilson loops, computed using Gaussian distribution for gauge fields (as is done by 
diagrams (d)-(g), for which only the two-point function of the gauge field enters), simply exponentiate: 
{e^ = exp(i(/ A J A)). As a consequence, the collision kernel as defined from (|9]) is tree-level exact 
in such theories; there is no interference between scattering events. 
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Part of diagram (f) is already included by the exponentiation of (ITU]) (diagram (a)): 
this generates the approximation to (f) in which the intermediate eikonal propagators 
are put on-shell. To avoid double-counting, this must be subtracted. We must first 
regulate the associated "pinching" (g^ — >■ 0) singularity, which we do by flowing a small 
external z-momentum uj into the Wilson loop. We then take the limit w — )■ after the 
subtraction. Diagram (g) poses no difficulty. 



(q-p) 



1 



{p)G++{q±-p) 



+ 



UJ — le 



{Pz 



(16) 



(17) 



Eg. ( !T6l) has a well-defined w — )■ limit, as follows from the identity l/{pz+ie — 
l/{pz—ie) = —27ii6{pz)- This limit takes a form identical to (ITTl) and the sum is 
proportional to Ca_, again as required by abelian exponentiation. This confirms the 
correctness of our evaluation of (f). 
The sum of diagrams (d)-(g) yields: 



C{qA 



){d)...(g) 



g^T^CsCA 



G+MG+^{qL-p) - 2G+MG++{qA_) 



{Pz 



le 



rriD 



47i{ql+ml) 



q^+Arri^) {q\+m 



(18) 
(19) 



The function G++ is Goo + Gzzi as given in eq. ([H]). To evaluate ( ITS]) we have found 
convenient to first apply integration by parts to the 1/ {pz — i^Y denominator, which re- 
moves the explicit p^-dependence and reduces the integral to a set of standard isotropic 
Feynman integrals. Eq. flTK]) is manifestly infrared- (and ultraviolet-) safe, upon enforc- 
ing p f-)- {q±—p) symmetry. 



4.5 Final formulae 

In summary, we have obtained all 0{g) contributions to the collision kernel C{qx): 
g'TCsml 



C(gx)(^^°) 
g'T^C^CA 



gi(gi+m' 



DJ 



7 



-uio — 2 



tan 



■niD 



32qi 



47r(g^+m|, 



2 \2 



tan-i ^ 

niu 



27cq±{ql + ml] 



rriD 



+ 



27rgi 



47r(gi+m2,) 



+ 



mn — '^^^^"^-p tan ^ 



DJ 



2m r 



(20) 
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These expressions are valid for ^ T. The leading order kernel for q± ^ T gets 
slightly modified, see fl22|) . 

The reader might wonder as to the appearance of arctangents with two distinct 
arguments in (l20l) . They can be understood by looking in the complex g^-plane: 
tan~^(g_i_/2mD) has a branch cut starting at qj_ = —Am'j-), exhibiting its origin from 
the exchange of a pair of two quanta of mass m/5 (longitudinal gluons). The branch cut 
of tan~^(g_|_/m£)), starting at q'j_ = —m'jy, arises from the exchange of one longitudinal 
and one transverse gluon. Both arctangents occur, since both of these pairs of states 
can be exchanged. Exchange of two massless quanta also occurs, and generates 1/ ^/q^- 
type of discontinuities instead of arctangents. 

5 Evaluation of g^^^o) 

The effective theory approach we have used so far is valid for q± <^ T. As mentioned 
in section 12. 3[ however, the momentum broadening coefficient q (second moment of 
C{q±)) receives contributions from all scales up to a process-dependent cut-off gmax- In 
this section we will assum 

Q'max T. 

To separate the soft and hard contributions to q, we find convenient to introduce 
an auxiliary scale q* obeying mn <^ q* '^T: 

The soft kernel C{q±Y°^^ is given by ([20]). The hard kernel C{qA_Y''"^ describes 
tree-level 2 — )■ 2 scattering processes against plasma constituents, with self-energy 
corrections omitted on the exchange gluon (since they represent only ~ g"^ corrections 
for q± ~ T). The large particle energy E ^ T guarantees that the Mandelstam 
invariants s ~ ET and —t = q]_ obey |t| ^ s, so that the relevant scattering matrix 
elements assume the universal (eikonal) form oc /t^ . The kinematics force q'^ = qz for 
the momentum transfer q. In fact, these processes are precisely described by the central 
cut of the (four-dimensional) diagram (b) of fig. |3l Performing the qz integration in 
the expression for the collision rate (as done in [17]; more details can be found in [26]). 
one obtains: 

Ciq±r'' = ^ / yf^^^ [2CAnB{p){l + nsip')) + ANfTjuMi^ - n^{p'))] , 

q± J i^^j p 

(22) 

^'^ In the context of jet quenching this is (parametricahy) equivalent to £■ T/g"^, with E the 
smallest energy of the participants. 
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with p,p' the initial and final momentum of the target particle; p' = p + qz, Qz = (f = 
^p^p^K regime q± <^ T , p' ^ p and fl22|) reduces (as it must) to the large q± 

limit of (dUD, C(gx) ~ ^^mf^C^T/gl. 

Integrating (122]) over g to obtain the hard contribution to fl21l) . and expanding it in 
powers of q* /T, yields: 



"hard 



NfTf 



+ 



67r 



log + II log (^) - 0.068854926766592 • • • + + • • • 

log (^) + log (^) - 0.072856349715786 ... + .. 



with the omitted terms being suppressed by (g*/T)^ or more. The quoted numbers are 
from [26] ; we have verified the first five significant digits by direct numerical integration 
of ( 122|) I^H The ~ g*/T term arises from soft bosons with p,p' <^T and can be obtained 
in the soft approximation nB{p),nB{p') — >■ T/p,T/p'; it is also given in [26]. 

The soft contribution to fl2T]) . e.g. the second moment of fl20]) . admits the expansion: 



^soft mlg^T , f q* \ g^T^CAmo 



n o log ^ + 

Cs 27r \iTiD J 27r 



q* 37r2 + 10-41og2 



16mo 167r 



(24) 



with the omitted terms being suppressed by powers of mo/q*- The q* dependence of 
( 123]) and IHM cancels out in their sum, as it must do, producing the advertised formula 
03]). This cancellation provides a rather nontrivial check on the calculation. 

The reader might inquire as to whether we have consistently included all 0{g) 
contributions to q. Taking q* ~ g^^'^T, for instance, the omitted terms ~ {q* /T)"^ in 
( 123]) might naively appear to be 0{g), suggesting contributions from other, omitted 
terms. 

Estimates of this kind can be misleading, however, because q* is not a physical scale 
in this problem. The matching region rriD <^ q* <^ T can be described equivalently 
using the low-energy description (EQCD) or the full theory, ensuring that q* always 
disappears from final expressions. This is seen explicitly for the leading truncation 
errors ~ q* /T in ( 123]) and (IM]) : instead of producing O^g^^"^) corrections, as one would 
naively expect setting q* ~ g^^^T, they cancel against each other and the leading 
correction is 0{g), not O^g^^"^). Since similar cancellations are bound to occur at all 
orders, this simply means that the scale q* should not enter power-counting estimates. 
Because higher loop diagrams are ~ g"^ when q± ^ T and because we have included 
all 0{g) effects when q± ~ rriD, we thus conclude that we have included all 0{g) 
contributions. 

Finally, we note that, in the spirit of j40j, we could have used dimensional regu- 
larization to separate the q integration, instead of the sharp cutoff q*. In this scheme. 



^^I thank P. Arnold for pointing to me a numerical error in an early draft of this paper. 
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the hard q* /T term in fl23p disappears: there is no suitable dimensionful parame- 
ter to replace q* . The 0{g) corrections then come solely from the (unambiguous) 
dimensionally-regulated soft integral ( 120|) . 



6 Jet Evolution 

We now extend the calculation to obtain the collision kernel relevant for bremsstrahlung 
and pair production processes. The new complication is that, except for QED processes, 
the relevant object to evolve in the plasma is no longer a "dipole": it involves three 
charges. For instance, to describe the gluon bremsstrahlung process ip — )■ gip, one must 
evolve an operator which annihilates a quark and creates a quark-gluon pair (see [2S] 
j29] |30] . which, however, use somewhat different notations): 

0^^^g = \^,g){^\. (25) 

The three color charges in fj25|) are paired together to form a color-singlet state, as 
dictated by the (DGLAP) gluon emission vertex which generates this operator. 

It turns out that only one transverse momentum suffices to describe the internal 
state of ( 125|) . A priori the description of three charges might seem to require three 
momenta, but momentum conservation reduces this by one, and a symmetry removes 
yet another one: by suitably choosing the z-axis it is always possible to "gauge" to 
zero one of the transverse momenta (see the discussion preceding eq.(6.6) in |3oj^ . 
In the following, for concreteness, we shall gauge to zero the transverse momentum of 
particle 1, and q±_ will refer to the transverse momentum of particle 2. 

At the leading order, the relevant collision kernel is a sum over two-body contribu- 
tions EHI |29] EOl: 



'^r3(g_L) C'2+C3-Ci ~ C'i+C3-C2 ~ , Ci+C2-C^ ~ El 
d-qJClnf = 2 + 2 ^^E^^^^ + 2 ^^E^^^^ ^'^^ 

with Ci and Ei respectively the Casimir and longitudinal momenta of the participating 
particles; C{q±) = ■^C{q±) denotes a single-particle collision kernel with its Casimir 
factor stripped off; we recall that the LO (and NLO) kernels respect Casimir scaling. 
In the special limit in which one of the Ei becomes much smaller than the other ones, 
the motion of this particle dominates and the kernel ( 126|) reduces to the one for single- 
particle diffusion, C{q±), for z = 2,3, and C{^q±) when i = 1. 

As we presently show, it turns out that the formula fl26l) also holds at NLO, provided 
the NLO expression ( l20i) for C{q±) is used in it. 



For high-energy jets (when at least one of the energy of the participant is large, E'max ^ T), 
these rotations can be taken to have energy-suppressed angles ~ q±/E, and thus to have negligible 
effects on the longitudinal momenta. Even when E^i^^ ~ T, the angles are at most ~ g and the 
changes in longitudinal momenta are ~ g^, beyond the accuracy considered in this paper. 
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(h) (j) 



Figure 4: Additional diagrams for the evolution of a triplet of charges. 
6.1 "Three-pole" propagation at NLO 

To keep the discussion simple we will assume that particle 3 is a gluon (color adjoint 
state), which is sufficient to cover all splitting processes in QCD (and JV = 4 super 
Yang-Mills). This ensures that particles 1 and 2 are antiparticles to each other. We 
denote by |s) the relevant singlet state in the tensor product of the three charges; 
explicitly, |s) is given by the representation matrices {ti)fj. 

The previously treated dipole diagrams (a)-(g) must now be summed over the three 
possible pairs of particles, and we must recompute their group theory factors. Diagrams 
(a)-(b) involve, in the case the interaction is between particles 1 and 2 |30j . 

which reproduces the structure (l26!l . upon summing over pairs and using rotational 
invariance to "gauge" to zero particle I's ±- momentum. Diagrams (c)-(g) fit the same 
structure, as follows from the fact that they organize themselves into commutators. 
For instance, 

(c) (X ^r»'^{s\tlt1®f,\s) = -^{s\t1®e,\s), 
if) + ig) (X {s\[t1,i]^^^\s) = -^{s\t-,0t-,\s). (28) 

Here we have used the identities [f'^t^] = ifH" and /"^c/'^^'^' = CaS""' . 

There are also new diagrams (fig. H]), which couple together the three particles 
nontrivially. In diagram (h) (see fig. S]), the Yang-Mills 3- vertex generates a factor /"^^ 
and the coupling to the gluon line is given by (t^)'^^ oc /^'^^, whence: 

(h) = {s\t1tltl\s)r'" oc Tr 1 {ftHH') p^^f-'^ = , (29) 

with the trace taken in the representation of the particle 1. We could prove this identity 
by making extensive use of the antisymmetry of the f°''"^. Diagrams (i) are similar to 
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diagram (g) treated in subsection \4A\ and the main point is that there is a sign between 
the two diagrams, due to the reversed middle propagator, thus yielding zero: 

(z) oc (s|t^ O [t^,t^] = 0. (30) 

Thus the new diagrams (h)-(i) vanish and the factorization formula f l26|) remains valid 
at NLO. 

We view this as somewhat surprising; this could be an artefact of the relatively low 
order in perturbation theory to which we are working. 

6.2 Discussion of operator ordering issues 

We now briefly discuss operator ordering issues, for the Wilson lines in iQ and their 
three-particle generalization ( 125|) . Although this is not directly relevant to the purely 
classical effects which are the main object of this paper, since nonperturbative defini- 
tions of q have been used in the literature [21] [19] we feel that a discussion of them 
can be of interest. 

To help clarify the physical significance of these issues, let us first consider, in QED, 
the processes of photon bremsstrahlung from a charge and of pair production from a 
photon. These processes differ in that the former takes place within the electromag- 
netic field generated by the initial charge, but the latter takes place in an essentially 
undisturbed medium (the induced field being suppressed by the small size of the pro- 
duced dipole). The collision kernels relevant to these two processes could thus be 
different, due to the different backgrounds, and should be defined differently. In the 
eikonal regime, it is the role of the Wilson lines trailing behind the charges to source 
the backgrounds, which requires that they be properly ordered. 

The proper ordering can be readily described using the language of the Schwinger- 
Keldysh "doubled fields" [H], in which amplitudes and their complex conjugate are 
described by type-1 and type-2 fields, respectively. For photon bremsstrahlung, evolv- 
ing the relevant matrix element requires one type-1 (and 7) and one type-2 
ijj field, whereas for pair production, evolving requires both charged fields to be 
type-1 (and 7 to be type-2). In the latter case the Wilson lines nearly cancel against 
each other (for a small dipole), whereas in the former case they fail to cancel, due 
to operator ordering issues (they live on different branches of the Keldysh contour): 
instead they source an electromagnetic field. This reproduces the expected physics. 

The story for QCD must be similar: for instance, evolving a \ip,g){il)\ operator, 
relevant for gluon bremsstrahlung, should require type-1 ip and g fields, and a type-2 
ip field, with the obvious replacements to be made for other processes. Thus we see 
that the strong coupling calculations of the momentum broadening coefficient in [21] 
and [22J, strictly speaking, gives a q applicable to photon bremsstrahlung, whereas 
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the "jet quenching parameter" defined in [T^], being defined from a space-hke hmit of 
correlators, is by hypothesis independent of operator ordering. 

It is not clear, at least to the author, the extend to which these effects can be 
numerically important. Obviously, at weak coupling, they are suppressed by a power 
of the coupling (the preceding subsection shows that the suppression is at least ~ g^). 
Furthermore, in the v 1 limit relevant to high-energy jets, an argument based on 
the shrinking down of the "causal diamond" enclosing any two points on the trajectory 
of the jet might suggest that these effects disappear — e.g. there is no time available 
for the induced field to influence the jet back again; a rigorous analysis, in particular 
of quantum effects, will not be attempted here. 
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A Relation to sum rule approach 

In this Appendix, we consider the problem of calculating, directly in four dimensions, 
the leading order collision kernel (|9]): 



with q± <^ T, and the full HTL-resummed propagator [11]. The simple result ( ITOl) 
for it has been obtained previously using a sum rule by Aurenche, Gelis and Zaraket 
(AGZ) [18], and our aim here is to establish the equivalence between our approaches. 

A.l Causality and Sum Rules 

Due to causality, retarded correlators G^{Q) must be analytic functions of the four- 
momentum Q when positive time-like or light-like imaginary four-vectors are added to 
it. This statement extends, in a Lorentz-covariant way, the familiar analyticity of 
in the g° upper-half plane [36] . 

Light-like imaginary parts are also allowed, because causality is preserved along 
light-fronts (e.g. G^{x~^) vanishes for negative light-cone time x^). 

In the classical approximation UB^q^) ~ T/q^, ( 1311) becomes: 




(31) 



(EH) = T 



/ 



dq.Gl^{q' = q.,q^)-Gi^{q^ 



(32) 



2vr q^ 
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To evaluate this by contour integration, we first move the Qz = pole slightly off- 
axis, l/g^ — 7- — ze), which does not change the result since (G^—G^) vanishes at 
Qz = 0. Next, we note, using the standard HTL expressions |T3], that vanishes 
(like I/Qz) at large |g^|, making it possible to close integration contours at infinity. 
Closing the contour for G^ (resp. G^) in the upper (resp. lower) half-plane, one 
obtains a unique residue iTG^^{q'^ = Qz = 0,q±) from G^ and nothing from G^, due 
to their aforementioned analyticity properties, thus reproducing ( ITOll : 

m = T(^- -T^] ■ (33) 

Additional poles at the Matsubara frequencies q^ = qz = 27rmT would have ap- 
peared in this result, in agreement with the sum (^^, had we kept the full Bose dis- 
tribution function nB{q^)- This shows that the classical approximation to distribution 
functions is equivalent to keeping only the n = Matsubara frequency. 

More generally, at higher orders in perturbation theory one could also imagine 
computing the dipole amplitude ([9]) using the real-time formalism. The fact that the 
operator of interest is supported at a constant value of x+ allows the integral over the 
conjugate momentum g"*" to be performed by contour integration, in the same way that 
integrals over g° (or equivalently, the sum-integrals of the imaginary-time formalism) 
would be done for equal-time correlators [36]. In addition to standard gauge field 
propagators, the integrands to be met also involve eikonal propagators ~ but 
since they do not depend on q^ they would not interfere with this calculation. Thus, as 
claimed in the text, we can also derive all the results of section H] using causality-based 
sum rules. This is in fact how we first obtained them. It is also clear that, in general, 
making the approximation nsip'^) — ?■ p^/T in all propagators in real-time is equivalent 
to dropping all n 7^ Matsubara modes. 

A. 2 AGZ's sum rule 

AGZ [18j study exactly the integral (!3T|) . but parametrized using a different variable, 
X = (f jq (so g'^(x) = qz{x) = \q±\x/V^ — x"^)'- 

dx 

m = \<lA J ^ 2vr(l-x^)3/2 ^++(^'^^) ^34) 

A key observation in [18] is that the HTL propagators, viewed as a function of x with 
g_L fixed and q^ = qz, are analytic in the whole complex x-plane, apart from a branch 
cut at real x G [—1, 1]. Using methods of complex analysis, they could then derive the 
result das]). 

To show that this analyticity property in x is equivalent to the analyticity in q'^ 
that we have used above (e.g. to causality), we rewrite the change of variable above 
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dMD as: 

= qz{x) = i\q±\^==, (35) 
Va; — 1 

and choose to put the branch cut of the square root at real x G [—1, 1]. Thus g° — t- i\q±\ 
as |x| ^ oo in any direction. This choice of branch cut ensures that G^{q'^{x), qz{x), q±) 
goes into the standard retarded function as Imx — )■ O"*", and is consistent with the 
conventions of [18], e.g. this function has the same analytic structure as the G^{x, q±) 
of [IE]. Careful inspection of ( l35l) then reveals that the imaginary part of g° is positive 
for all X, establishing that analyticity in x (for q^ = q^ and at fixed q±) is a consequence 
of the above-discussed analyticity in g+. It thus applies to any propagator, extending 
the claim of jlSj . 

The authors of [18] worked in the Coulomb gauge and found, at intermediate steps, 
contributions from the large circle at |x| = oo (proportional to l/(g^+|m|,)), which 
in the end, precisely canceled out between the longitudinal and transverse channels. 
Since no such term has appeared in our approach, the reader might wonder as to 
their claimed equivalence. What happens is that these contributions are mere gauge 
artefacts; this also explains their ultimate cancellation. To see this, we note that the 
residue at |x| — )■ oo corresponds to a pole at q^ = qz = iq±, which is at an ordinary 
point in the upper-half g^-plane and is thus forbidden by causality. But since a gauge 
like the Coulomb gauge does not respect causality in a Lorentz-covariant sense (its 
field mediates an instantaneous Coulomb interaction), such poles are not forbidden in 
individual, gauge-dependent terms. They are bound, however, to cancel out in physical 
quantities like C{q±). Our approach assumes Lorentz-covariant causality from the start 
and cannot detect such unphysical contributions. 

These complications, associated with causality violations in non-covariant gauges 
(in intermediate expressions), are easily avoided by working in a covariant gauge. 
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